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Abstract 
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timal performance. Recent regulations demand the use of environmentally un¬ 
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1 Introduction 


Many currently used refrigerant gases have a great impact on the environment since 
they either contribute largely to the greenhouse gas effect, or because they tear the 
ozone layer, or both. In an attempt to protect the environment, regulations preventing 
the production and use of certain refrigerant gases have been implemented ||T|. 

Gas detectors are widespread for detection, tracking and triggering of charged parti¬ 
cles such as muons in Nuclear and High Energy Physics (HEP). They are characterized 
by simple and reliable use, but utmost care must be taken to issues such as properties 
of gas interaction with materials, gas purification, gas mixture contaminants, etc. O- 

m. 

A large part of gas muon detectors used in HEP operates with mixtures containing 
the regulated refrigerants as quenching medium in applications where excellent time 
resolution and avalanche operation are necessary. Therefore, actions towards finding 
new mixtures must be undertaken. Gas Electron Multiplier (GEM) m detectors oper¬ 
ate in experiments such as CMS (Compact Muon Solenoid) at the LHC (Large Hadron 
Collider) with an argon/C02 mixture urn . However, for high time resolution applica¬ 
tions an argon/C 02 /CE 4 mixture is used lIT^ . where CE 4 has a Global-Warming Poten¬ 
tial (GWP) of 7390 IIT 3 II . Resistive Plate Counters (RPC) ifTOl currently operate with 
a E-based R134a/Isobutane/SE6 gas mixture, with typical GWP of 1430. Investiga¬ 
tions into new gas mixtures has to be performed in order to keep the mixture properties 
while complying with the regulations. A few industrial refrigerant industrial replace¬ 
ments were proposed ifTHl as alternatives to R134a. A study of transport properties of 
currently used gas mixtures in HEP, and evaluation of transport properties of freon-less 
gas mixtures, was recently published ITSlllfT^ . 

The aim of this paper is to discuss some of the important properties of gases for 
particle gas detectors, to list and summarize basic properties of eco-friendly refriger¬ 
ants from the literature, to discuss their properties for materials compatibility and safe 
use, and make a prediction on selected parameters crucial for the performance of gas 
detectors considered. While this study is aimed to GEM and RPC detectors, its findings 
can be considered for selection of ecogas replacement for other gas detectors. 

2 Gas properties 

Eor a gas mixture to be appropriate in a gas detector, first of all it has to comply with 
the regulations. Eurthermore, its properties must also be appropriate for the specific 
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type of detectors. For example, a gas that is suitable for the RPC detectors may not be 
fully optimized for the GEM detectors. To better find the appropriate gas for a detector, 
an understanding of the infiuence of different parameters is required. This section aims 
to clarify the most essential parameters for gases. 

2.1 Global-warming potential and Ozone depletion potential 

In order to estimate the impact of a refrigerant on the environment, the effects have 
to be quantified. Two important effects are the contribution to the greenhouse effect 
and the depletion of the ozone layer. The first mentioned effect is measured in Global- 
Warming Potential (GWP), and is normalized to the effect of CO 2 (GWP =1), while the 
effect on the ozone layer is measured in Ozone Depletion Potential (ODP), normalized 
to the effect of CCI 3 F (ODP =1). The effects of selected refrigerant candidates are 
listed in table [T] 

2.2 Stopping power 

When a particle passes through a medium, energy is transferred from the particle to 

the surroundings. The energy lost is typically defined as the stopping power expressed 

as - ( —— ), where p denotes the density of the medium, E denotes energy, and x 
p \ ax / 

is length. This expression is itself dependent on the momentum p of the particle. In 
the next section it will be shown the particle will experience a minimum energy loss 
at around pjMc ^ 3, where M is the particle mass, and c is the speed of light. By 
knowing the average minimum energy loss, we also know the minimum energy the 
particle loses by passing through a detector. The minimum mean ionization energies 
for the refrigerants under consideration are summarized in table 

2.3 Radiation length 

The radiation length Xq is a characteristic length of a medium. It describes both the 
mean distance required for a high energy electron to lose all but e~^ of its energy 
due to bremsstrahlung, and of the mean free path of a e+e“ pair produced by a 
high-energy photon ca. These quantities are also estimated and summarized in table 

El 

2.4 Ion pair production 

When an incoming particle passes through a medium, it will eventually interact with 
the medium and transfer some of its energy to ionize atoms. In this process, a pair con- 
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sisting of an ionized atom and a free electron is produced. The number of ionizations 
produced by an incoming particle per unit length is denoted by Np, in units of cm“^. 
Each produced ion pair will have an initial kinetic energy and can itself produce an ion 
pair, called secondary ion pair production. The sum of the primary and secondary ion 
pairs production per unit length is denoted Np, and will be mainly depending on the 
material and the incoming particle energy and mass. This parameter is relevant in parti¬ 
cle gas detectors as it determines both the number and the size of avalanches produced 
by a single incoming particle when the gas is under an amplifying electric field. 

2.5 Drift velocity 

In the absence of electric field, electrons move randomly in all direction having an av¬ 
erage thermal energy 3/2KT. In presence of an electric field, the electrons start to 
drift along the field direction with mean drift velocity Vd (the average distance covered 
by the drift electron per unit time). The drift velocity depends on the pressure, temper¬ 
ature and could be affected by the presence of pollutants like water or oxygen. Electric 
field lines are, in addition, deformed by any magnetic field applied to the gas volume. 
Einally, electronegative pollutants deplete the gas of electrons. Computer codes such 
as Magboltz CD computes the drift velocity using montecarlo techniques, by tracing 
electrons at the microscopic level through numerous collisions with gas molecules. 

2.6 Diffusion 

Electrons and ions in a gas are subject only to an electric field and move on average 
along the electric field, but individual electrons deviate from the average due to scatter¬ 
ing on the atoms of the gas. Scattering leads to variations in velocity, called longitudi¬ 
nal diffusion, and to lateral displacements, called transverse diffusion. The scattering 
process in each direction can, to a good approximation, be considered gaussian on a mi¬ 
croscopic scale. In cold gases like carbon-dioxide for example, the diffusion is small, 
while the drift velocity is low and unsaturated for values of electric fields common in 
gas detectors; this implies a non linear space time relation. Warm gases, like argon for 
instance, have a higher diffusion; when mixed with polyatomic/organic gases having 
vibrational thresholds between 0.1 and 0.5 eV, diffusion is reduced in most cases, while 
the drift velocity is increased. 
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2.7 Townsend coefficients 


The average distance an electron travels between ionizing collisions is called mean 
free path and its inverse is the number of ionizing collision per centimeter a (the first 
Townsend coefficient). This parameter determines the gas gain of the gas. If no is the 
number of primary electron without amplification in uniform electric field, and n is 
the number of electrons after distance x under avalanche condition, then n is given by 
n = noe^^ and the gas gain G is given by G = no/n = The first Townsend 
coefficient depends on the nature of the gas, the electric field and pressure. To take 
into account the augmented emission of electrons by the cathode caused by impact of 
positive ions, it is customary to introduce r], Townsend’s second ionization coefficient 
or attachment parameter, i.e., the average number of electrons released from a surface 
by an incident positive ion, according to the formula 

^ax 

1 - - 1 ) 

2.8 Lorentz angle 

Because of the defiection effect exerted by a magnetic field perpendicular to the electric 
field and the motion of the electron, the electron moves in a helical trajectory resulting 
in a lowered drift velocity and transverse dispersion. Thus the arrival time of electrons 
at the anode changes and the spread in the drift time increases. The angle which the 
drifting electron swarm makes with the electric field is defined as the Lorentz angle 
of the particular gas or gas mixture under consideration. The Lorentz angle depends 
on both the electric field and the magnetic field. It is normally large at small electric 
fields but falls to smaller values for larger electric fields and is approximately linear 
with increasing magnetic field. 

2.9 Hazards and flammability 

Many refrigerants may constitute danger for the user and its environment. The great¬ 
est dangers involved are the fiammability and toxicity. We have used two standards 
in categorizing the refrigerants in table and The Ashrae standard ifT^ gives each 
refrigerant a number denoting fiammability from 1 (not fiammable) to 3 (highly fiam- 
able), as well as a letter A (non-toxic) or B (Toxic). The Health Material Hazardous 
Material Information System (HMIS), rates Health/ Flammability/ and Physical haz¬ 
ards from 0 (low) to 4 (high). 


5 



2.10 Compatibility with materials 

Some refrigerants are incompatible with certain materials, and can either react vio¬ 
lently, or have long term effect. Some refrigerants may even produce toxic decompo¬ 
sition and/or polymerization. Known incompatibilities and toxic byproducts are sum¬ 
marized in Table |3] and m 

2.11 Aging and longevity 

Aging is defined (following |[20l) the general deterioration of the detectors during their 
operation. The aging phenomenon is very complex and depends on several parame¬ 
ters. The commonly used variables include the cross-sections, electron/photon ener¬ 
gies, electrostatic forces, dipole moments, chemical reactivity of atoms and molecules, 
etc. For a comprehensive (although non recent) collection see EDEaCa. A more 
recent review of ageing effects in GEM detectors can be found in 1^ . 


3 Estimation of Gas Parameters 

3.1 Stopping power 

Quantities such as the minimum ionization energy can be found if the stopping power is 
known. An approximate expression for moderately relativistic particles in the momen¬ 
tum region 0.1 < ^^7 = pjMc < 1000 can be found using the Bethe-Bloch equation, 
given by lf30l 


p\ dx / A 


dE 


1 , 2meC^P^YT, 


2 ^" 


J2 


-/? - 


m) 


( 2 ) 


where ( —— ) is the mean energy loss per length, p is the density of the medium, 
\ dx / 

I is the mean excitation energy, and is the density effect correction function 

to ionization energy loss. iT is a constant given by AirNArlrrieC^, and Tmax is the 
maximum energy transfer in a single collision, given by 

2meC^/5^7^ 


T — 

max — 


(3) 


1 + 2^melM + {melM)‘^ ’ 

where M is the mass of the incoming particle. 

The mean excitation energy / for a composite medium can be approximated from the 
composite atoms by the relation ISTIl . 


/ = exp 


rcj {ZjIAj) In Ij 


/(^M) 


(4) 
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Molecular name 

Chemical 

formula 

CAS 

Refrigerant 

identifler 

GWP 

ODP 

Chloropentafluoroethane 

C 2 CIF 5 

76-15-3 

R115 

7370 

0.44 

Hexafluoroethane 

C 2 FQ 

76-16-4 

R116 

- 

- 

2,2-Dicloro-l,l,l- 

trifluoroethane 

C 2 HCI 2 F 3 

306-83-2 

R123 

120 

0 

l-Chloro-3,3,3- 

Trifluoropropene 

C 3 H 2 CIF 3 

2730-43-0 

R1233zd 

4.7-7 

0 

2,3,3,3 -Tetrafluoropropene 

C 3 H 2 F 4 

754-12-1 

R1234yf 

4 

0 

1,3,3,3 Tetrafluoropropene 

C 3 H 2 F 4 

29118-24- 

9 

R1234ze 

6 

0 

Trifluoroiodomethane 

CF 3 I 

2314-97-8 

R13I1 

0 

0 

1,1,1,2-Tetrafluoroethane 

CH 2 FCF 3 

811-97-2 

R134a 

1430 

0 

Tetrafluoromethane 

CF 4 

75-73-0 

R14 

7390 

0 

1,1,1 -trifluoroethane 

CH 3 CF 3 

420-46-2 

R143a 

4300 

- 

1,1 -Difluoroethane 

C 2 H 4 F 2 

75-37-6 

R152a 

124 

0 

Octafluoropropane 

C 3 F 3 

76-19-7 

R218 

- 

- 

Propane 

C 3 H 8 

74-98-6 

R290 

3 

0 

Difluoromethane 

CH 2 F 2 

75-10-5 

R32 

650 

0 

Isobutane 

C 4 H 10 

75-28-5 

R600a 

3 

0 

Sulfur Hexafluoride 

SFq 

2551-62-4 

R7146 

23000 

0.04 

Carbon Dioxide 

CO 2 

124-38-9 

R744 

1 

0 

Octafluorocyclobutane 

C 4 F 8 

115-25-3 

RC318 

- 

- 

Pentafluoroethane 

HF 2 CF 3 

354-33-6 

R125 

3400 

0 

Trifluoromethane 

CHF 3 

75-46-7 

R23 

0 

0 

R409 : 

CHCIF 2 

75-45-6 

2837-89-0 

75-68-3 

R22 

(60%), 

R142b 

(25%), 

R124 

(15%) 

1700 - 620 

0.5 / 0.065 
/0.02 

R407c : 

CH 2 F 2 , 

CF 3 CHF 2 , 

CH 2 FCF 3 

75-10-5, 

354-33-6, 

811-97-2 

R32 (21- 

25%), 

R125 

(23-27%), 

R134a 

(50-54%) 

650 3400 

1430 

000 


Table 1: Summary of various refrigerant candidates. 
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Name 

I 

■(-) 

/min 

^0 

Np 


[eV] 

MeV-^ 

cm"'- 


r g 1 

-cm^- 

[cm“^] 


R32 

89.3602 

1.80973 

35.4581 

49.2 

R7146 

127.401 

1.67833 

28.6027 

92.0 

R600a 

47.848 

2.24057 

45.2260 

81.0 

R1234yf 

91.9674 

1.7734 

35.8204 

89.5 

R152a 

78.1889 

1.88706 

37.0969 

67.1 

R1234ze 

91.9674 

1.7734 

35.8204 

89.5 

R115 

116.695 

1.69178 

29.2197 

98.4 

R1233zd 

106.689 

1.73915 

29.7636 

105 

R290 

47.0151 

2.26184 

45.3725 

65.2 

R13|l 

271.737 

1.42486 

11.5399 

272 

R134a 

95.0294 

1.76439 

35.1542 

81.6 

R14 

107.127 

1.69909 

33.9905 

63.6 

R123 

125.275 

1.69722 

25.5416 

98.4 

R143a 

87.8152 

1.8126 

35.8928 

74.8 

R744 

88.7429 

1.81124 

36.1954 

37.2 

R23 

99.9508 

1.7402 

34.5214 

56.9 

R116 

105.075 

1.70566 

34.2947 

93.3 

RC318 

101.578 

1.71721 

34.8435 

123 

R218 

104.13 

1.70873 

34.439 

117 


Table 2: Minimum ionization, radiation length and number of primary ion pair creation 
for the considered refrigerants, as well as the approximated mean ionization energy 
used. 



Refrigerant 

Molecular 

Density g/L 

Boiling point 

HMIS 

Ashrae 


weight 


degC 


Safety 

Group 


Material incompatibility 


Hazardous decomposition products and polymerization 


R115 


154.4 

6.623 

-39.1 

1/0/2 


A1 


Material is stable. However, avoid open flames and high temperatures. Incompatible with alkali or 
alkaline earth metals-powdered Al, Zn, Be, etc. 


Decomposition product are hazardous .’’FREON” 115 Fluorocarbon can be decomposed by high tem¬ 
peratures (open flames, glowing metal surfaces, etc.)forming hydrochloric and hydrofluoric acids, and 
possibly carbonyl halides. Thermal decomposition can yield toxic fumes of fluorides such as Hydrogen 
Fluoride, Hydrogen Chlo- ride. Carbon Monoxide and Chlorine. 


R116 


138.01 

5.734 

-79 

1/0/0 


Al 


May react violently with alkaline-earth and alkali metals Thermal decomposition yields toxic products 
which can be corrosive in the presence of moisture. If involved in a Are the following toxic and/or 
corrosive fumes may be produced by thermal decomposition: Carbonyl fluoride. Hydrogen fluoride. 
Carbon monoxide 


Thermal decomposition products: halogenated compounds, oxides of carbon 


130.5 

6.10 

-19 

- 


R1233zd 


Incompatible with polyacrylate, Viton®, natural rubber, silicon rubber and other elastomers. 


Is considered non-toxic at less than 800 ppm[? ]. Hazardous polymerization can occur. If involved in a 
fire, production under thermal decompose into pyrolysis products containing Hydrogen Fluoride, Carbon 
Monoxide, Carbonyl halides, and Hydrogen Chloride can occur E). 


R1234yf 


114.0 

4.82 

-29 

0/2/2 


Incompatible with alkali metals, Zn, Mg and other light metals. 


If involved in a fire, production under thermal decompose into pyrolysis products containing Fluorine, 
Carbon Monoxide, Carbonyl halides, and Hydrogen halides can occur. No toxic decomposition should 
happen under normal conditions. 


R1234ze 


114.0 

4.82 

-29 

1/0/0 


Incompatible with strongly oxidizing materials and finely divided Mg and Al. 


If involved in a fire, production under thermal decompose into pyrolysis products containing Fluorine, 
Carbon Monoxide, Carbonyl halides, and Hydrogen halides can occur. Polymerization may also occur. 


R13I1 


195.9 

- 

-22.5 

- 


Incompatible with active metals, fires of hydrides, and materials containing oxygen. 


Can decompose to Iodine, Hydrogen Fluoride, and Hydrogen Iodide. 


R134a 


102.0 

4.320(501 

-26.5 

1/1/0 


Al 


Chemically reactive with K, Ca, powdered Al, Mg, Zn. Under high temperature/ high pressure, may 
react with Al surfaces. 


Under special circumstances (e.g. high temperature) Carbon monoxide. Carbonyl fluoride. Hydrogen 
fluoride can be produced. Under normal storage and use, no hazardous decomposition should be pro¬ 
duced 


88.0 

3.65 

-128 

0/0/0 


R14 


Al 


Not compatible with aluminium, alloys containing more than 2% magnesium, alkali metals in powdered 
form and carbon dioxide above 1000 °C. 


If involved in a fire, production under thermal decompose into pyrolysis products containing hydrogen 
fluorid and carbonyl fluoride. 


Table 3: Chemical, physical and compatibility information of the refrigerants (Part 1). 
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Refrigerant 

Molecular 

Density g/L 

Boiling point 

HMIS 

Ashare 


weight 


degC 


Safety 

Group 


Material incompatibility 

Hazardous decomposition products and polymerization 


R143a _ I 84.0 _^_ | -47.6 __^2_ 

Can form explosive mixture with air. May react violently with oxidants. Air, Oxidiser. Non recom¬ 
mended: Hydrocarbon based lubricant, significant loss of mass by extraction or chemical reaction and 
Fluorocarbon based lubricant, significant loss of mass by extraction or chemical reaction. 

Thermal decomposition yields toxic products which can be corrosive in the presence of moisture. 

R125 I 120 I 1.24g/cm3 | -48.5 | 1/1/0 | A1 

Under very high temperature and/or appropriate pressures freshly abraded aluminum surfaces may cause 
strongly exothermic reaction. Chemical reactive metals: potassium calcium powdered aluminum, mag¬ 
nesium and zinc. 

The product is stable. Do not mix with oxygen or air above atmospheric pressure. Any source of high 
temperatures, such as lighted cigarettes, flames, hot spots or welding may yield toxic and/or corrosive 
decomposition products. 

R22 I 86.45 | 3 | -40.1 | 1/0/0 | - 

Chemically reactive metals: potassium, calcium, powdered aluminum, magnesium, and zinc, powdered 
metals, powdered metal salts 

Hazardous decomposition products: Halogens, halogen acids and possibly carbonyl halides. Carbon 
monoxide. Phosgene, Hydrogen chloride. Hydrogen fluoride. Carbonyl fluoride. 

R744 I 44 I 1.52 | -78.5 | 1/0/0 | A1 

The product is stable under regular conditions 

Materials to avoid: strong oxidising agents, strong acids. Hazardous decomposition products: In com¬ 
bustion emits toxic fumes. 

R142b I 100.5 I 4.18 | -10 | 2/4/0 | - 

Materials to avoid: Light and/or alkaline metals. Alkaline earth metals. Powdered metals. Oxidizing 
agents. Chlorine, Powdered aluminum, magnesium, zinc, beryllium and their alloys. 

Hazardous decomposition products: Gaseous hydrogen fluoride (HF)., Gaseous hydrogen chloride 
(HCl)., Fluorophosgene, Phosgene 

R152a I 66.1 | 2.738 | -25 | 1/4/2 | A2 

Extremely reactive with oxiding materials , such as alkaline, alkaline earth metals, and other reactive 
chemicals, (i.e. Na, K, Ca, Mg, powdered Al, Zn), brass, and steel . Incompatible with amines, bases, 
and halogens. 

Under normal condition, hazardous decomposition and/or polymerization products should not be pro¬ 
duced. If exposed to fire, hazardous products may be produced. 

R218 I 188.0 I 8.17 g/1 gas | -36.7 | - | Al 

Stable under normal conditions materials with which gas mixture is incompatible: oxidizing materials 
and alkali and alkali earth metals. May react violently with chemical active metals as sodium, potassium 
and barium powdered magnesium, powdered aluminum and organometallics 

Thermal decomposition yields toxical products which can be corrosive in presence of moisture hazardous 
decomposition products: acid halides 

R23 I 70.0 I -2.946 kg/m3 | -84.4 | - | Al 

Incompatible materials: metals,poly styrene, natural rubber, alloys of more than 2% magnesium in the 
presence of water, nitrosyl fluoride,N 203 , lime at dull red heat, and metals at elevated temperature 

Decomposition products: halogenated compounds, oxides of carbon, hydrogen fluoride, thermal de¬ 
composition may produce toxic fumes of fluorides. Decomposition products may include the following 
materials: carbon dioxide carbon monoxide halogenated compounds carbonyl halides. 

Table 4: Chemical, physical and compatibility information of the refrigerants (Part 2). 
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Refrigerant 

Molecular 

Density g/L 

Boiling point 

HMIS 

Ashare 


weight 


degC 


Safety 

Group 


Material incompatibility 


Hazardous decomposition products and polymerization 


44.1 

1.86 

-42 

1/4/0 


R290 


A3 


Incompatible with acids, oxygen, oxidizing materials, copper, some plastics. Chlorine Dioxide. 


Under normal condition, hazardous decomposition and/or polymerization products should not be pro¬ 
duced . May produce carbon monoxide and other toxic gasses under thermal decomposition. 


52.0 

11.4 

-51.7 

1/4/1 


R32 


A2 


incompatible with acids and oxidizing materials as Na, K, Ca, Zn, Mg, powdered Al, and other active 
metals. Incompatible with air and moisture. 


No hazardous decomposition/polymerization should be produced under normal conditions. 


58.1 

8.93 

-11.7 

1/4/0 


R600a 


A3 


Incompatible with oxiding materials, halogenated hydrocarbons, halogens, and metal catalysts. 


No hazardous decomposition/polymerization should be produced under normal conditions. May produce 
carbon monoxide and other toxic gasses under thermal decomposition. 


146.1 

6.17 

-63.7 

1/0/0 


R7146 


Stable with most chemical, except metals other than aluminium, stainless steel, copper brass, silver, at 
elevated temperatures (^204° C). Also reacts violently with disilane. 


Decomposes into Sulfur oxides and hydrogen fluorine. 


Table 5: Chemical, physical and compatibility information of the refrigerants (Part 3). 


where Wj , Zj , Aj and Ij is the fraction by weight, atomic number, atomic weight and 
mean ionization energy, respectively, of the j’th constituent. The shape of the 5 {^5) 
function for non-conducting materials can be approximated by 

ifX > Xi; 

6{-f6) = \ 2(lnlO)X-C ifXo<X<Xi; (5) 

[ 0 ifX<Xo; 

where X = logio(7^)- finding an approximate expression for the parameters C, 
Xq and Xi based on experimental fits, we refer to 1521 . For gases with momenta below 
the density effect correlation function can be neglected. A plot of the calculated 
energy loss for different refrigerants is show in figure 


3.2 Radiation length 

The radiation length of an atom can be found by ca 

Vo = 716.405(cm-2mol)^/ [Z\Lrad - f{z)) + (6) 

f{z) = ( 2 > 2\ ~ 1-2022 - 1.03692^ + ——— (7) 

n=l ^ ' 

where Li and L 2 is given by table |3.2| f{z) is the one-photon exchange approxima¬ 
tion, and ^ = aZ, a being the fine-structure constant and Z is the atomic number. 
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Figure 1: Energy loss as a function of the relativistic time dilation factor for various 
refrigerants. 
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z 


Li 


L2 


1 

2 

3 

4 

5> 


5.31 

4.79 

4.74 

4.71 


6.144 

5.621 

5.805 

5.924 


5> ln(184.15Z-i/3) ln(1994Z-2/3) 

Table 6 : Li and I /2 expressions for various atom numbers. 


This formula, however, only holds for free atoms. When one wants to find the stopping 
power for a molecule, one has to take into account the influence from molecular bind¬ 
ings, crystal structures and polarization of the medium. By neglecting these effects, 
however, one can find an approximate expression by weighting the radiation length of 
the single atoms 



( 8 ) 


where j reffers to the j’th atom’s constituent. 

3.3 Estimation of ionization pair production 

In order to model the number of primary ionization caused by a single particle, the 
cross section for all the particle-atom interaction should be calculated. The number 
of primary electrons per unit length would then be the integral over energy across all 
the energy transfer cross sections. This is tedious work since all electron orbital trans¬ 
fers have to be considered. An easier, but approximate correlation between primary 
ionization and atom number has been found based on experimental data by (sa 



(9) 


which holds for normal pressure and temperature (NPT) (1 atm, 20°C). For different 
pressure and temperature, the number scales with the density. This value should only 
be taken as a rough estimation though. This formula has proven to work best for hy¬ 
drocarbons and worst for molecules consisting mainly of fiuorine, differing as much as 
30% from the experimental value for CF 4 . 

The total number of ionization has proven to be more difficult to estimate. Whereas 
no general formula has been derived, the most straightforward method will be to use the 
cross sections used to calculate the primary ionization electrons, and use Monte Carlo 
simulations to track the production of secondary electrons from primary electrons. The 
total number of pair ionization turns out to be dependent on the incoming particle 
energy and mass, and a general expression can therefore be difficult to find. For an 
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incoming particle, W = defines the average energy cost to produce an ion pair. 
Luckily, this turns out to be a function slowly varying with particle energy (271, and can 
therefore be taken to be a constant in an energy interval. The total ionization per unit 
length can then be found by 



( 10 ) 


If the W values for specific gases are know, the average W value for a gas mixture can 
be found by (531 



( 11 ) 


where rim denotes the index of the molecule, and /(n^) denotes the relative number 
of molecules of the given sort in the mixture. 

The value of W is difficult to predict, and there is not a direct way to give a proper 
estimate based on experimental data alone. A montecarlo simulation is in preparation 
which uses the photoabsorption ionization and relaxation (PAIR) model (5^ and it will 
be the subject of a forthcoming paper. 


4 Quantum chemistry calculations 


In this section of the paper we describe the calculation of gas properties using the quan¬ 
tum chemistry calculation tool NWChem . The ionization energies, electron affinities 
and the absorption spectrums are calculated for several kinds of potential Eco-friendly 
gas replacements. The results are validated and compared with different methods or 
calculation schemes. In the last part, the problems related to the computation of the 
Townsend parameter and prop sects for future works are discussed. 

Developments of quantum chemistry in the past few decades have improved the 
calculation to be accurate enough to be compared with real experiments. These calcu¬ 
lations enable us to study gas properties from atom level, where the underlying theory 
is believed to be more fundamental and may therefore provide us more general descrip¬ 
tions of the chemical properties. 

In the paper, we used the NWChem software (54l for computation, and the Py- 
Mol (55l tool for visualization. The calculation determines the electronic structure 
of the molecule system and can provide basic information of the system like ground 
state energy, dipole momentum, charged density and molecular orbit. To obtain prop¬ 
erties more suitable to be compared to experimental data, the software also provides 
frameworks to compute the absorption spectrum and the ionization energy. 
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Here we first show the basic geometrical information for interested molecules. 
Then we describe the procedure to calculate the absorption spectrum and in the follow¬ 
ing we present the methods to calculate the ionization energy and the electron affinity. 
In the end, we discussed the possibility to estimate the Townsend parameter under some 
certain environment of particle detectors. 

4.1 Molecules and their optimized geometries 

The freon gases we have been interested in are R134a (CH 2 FCFS), R152a (C 2 H 4 F 2 ), 
HF01234ze (CFHCHCFs), HF01234yf (CH 2 CFCFS), CF^I and HF01233zd 
(CHClCFCFs). Meanwhile we choose Rl2{CCl2F2) to be the standard freon gas 
model. CH^ and CF 4 ^ are the molecules we were used to make a comparison be¬ 
tween the NWChem calculated results and experimental results, in order to check the 
reliability of NWChem . Fig.Q shows the optimized ground state geometries of gas 
molecules, where the green balls stand for carbon atoms, the grey balls stand for hydro¬ 
gen atoms, the indigo balls stand for fiuorine atoms, the brown balls stand for chlorine 
atoms, and purple ball for iodine. Fig.([^ shows the highest occupied molecular orbitals 
(HOMO) of gas molecules. 

5 Calculation of absorption spectrum 

The excitation energy of a molecule is one of the fundamental properties of molecular 
interactions one can get from experiment. To study these properties, we used the frame¬ 
work introduced in Ref. tSbl to simulate the time-dependent response of molecules un¬ 
der external fields using quantum chemical calculations. The framework implemented 
is the Real-Time Time-Dependent Density Functional Theory (RT-TDDFT) ||57]| method 
within the NWChem , making it capable of doing the simulation beyond small pertur¬ 
bation from the ground state. 

In our case, we are mostly interested in low excitations of small molecules. We 
therefore adopt the procedure as described in section 3 of the reference 15^ . The 
choice of external field is (5-function-like electric field “kick” 

E(t) = K, • exp [—{t — /2w^] (i, (12) 

where to is the center of pulse, w is the pulse width, which has dimensions of time, 
(i = X, i is the polarization of the pulse, and k. is the maximum field strength. 
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A 
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(a) CH 4 

(b) CFa 

(c) RU(CCl 2 F 2 ) 

(d) CF3I(CF37) 



(e) R134a(Cif2FCF3) (f) R152a(C2if4F2) 



(g) 

HFO1 234zq{C FHCHCFs ) 



(h) 

HFO 1234yf (Cif 2 CFCF 3 ) 



(i) 

HFO 1 233zd{C HCICFCF 3 ) 


Figure 2: Optimized ground state geometries of gas molecules. 
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(a) CHa 


(b) CF 4 


(c) K\ 2 {CCl 2 F 2 ) 


(d) CF3I(CF3/) 



(e) R134a(Cif2FCF3) (f) R152a(C2if4F2) 



(g) 

HFO1 234zq(CFHCHCF3 ) 



(h) 

HFO 1234yf (Cif 2 CFCF 3 ) 



(i) 

HFO 1 233zd{C HCICFCF 3 ) 


Figure 3: Highest Occupied Molecular Orbitals (HOMO) of gas molecules. 


The system is then evolved in time, and the dipole moment can be calculated with 
respect to the added dipole coupling term 


. E{t), (13) 

where D is the transition dipole tensor of the system. Then we Fourier transform the 
dipole signals to construct the complex polarizability tensor and finally the 

dipole absorption spectrum is 

'5(w) = iTr[cr(w)] = ^Tr[Im[aij{uj)]]. (14) 

We first validate our calculation by repeating the calculation of CH 4 lowest ex¬ 
citation energy in Table 0 of the paper 1561 . With the same basis set (6-3 IIG) and 
functional (B3LYP), our calculation gives 11.16 eV excitation energy, which is consis¬ 
tent with the 11.13 eV in the paper. The small difference can come from other minor 
uncertainty sources like the choice of time separation in the simulation. 

We further checked the dependence of wave function basis sets as well as density 
functional of our results. The dependence of basis sets is found to be much smaller 
than the functional one. Fig. 0 and Fig. 0 compares the results using two different 
functionals. As shown in the figures, although the strength varies large between the 
two functionals, the excitation energies are similar, especially the lowest one, which is 
also the same story as in the reference ll56l . 

The same calculation is performed to gas molecules of our interest. The optical 
absorption spectrums are compiled in Fig. 0 - 

In the Fig. 0 we put together all spectra of interesting molecules. Comparing to 
the R12, which is one of the Freon gas spectrum, we found the similarity between 
molecules are in good agreement with their structures. In other words, the absorp¬ 
tion spectrum is closer to that of R12 if the candidate molecule is similar to the R12 
molecular structure. 

6 Calculation of the vertical ionization energy and the 
electron affinity 

The vertical ionization energy and the electron affinity are another two important pa¬ 
rameters closely related to experiments. In the case of small molecules, there is a 
simple way to do the estimation. The Koopman’s theorem states that at in closed-shell 
Hartree-Fock (HF) theory, the negative value of the Highest Occupied Molecular Or¬ 
bital (HOMO) energy is a good approximation of to the experimental ionization energy 
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(d) X-dipole, B3LYP 


(e) Y-dipole, B3LYP 


(f) Z-dipole, B3LYP 


Figure 4: R12(CC/2^2) dipole moment in the 3 kicked directions, using the 6-31 G* 
basis and LDA (top) and B3LYP (bottom) functionals. 



Figure 5: A comparison of R12(CC/2^2) absorption spectrum using different func¬ 
tionals. 
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Figure 6 : Absorption spectrum of interested gas molecules. The calculation is per¬ 
formed in the AO basis, where the solid line corresponds to the overall dipole strength 
and the dashed lines correspond to polarizations in the coordinate basis directions. 
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Figure 7: A comparison of absorption spectrum among molecules. 


21 























(IE). Similarly, it also suggests that the electron affinity (EA) is equal to the negative of 
the Lowest Unoccupied Molecular Orbital (LUMO) energy. However, the EA part of 
the prediction is in general not reliable because of the large effect of orbital relaxation 
on the LUMO eigenvalue 15^1591 . 

In this section, besides estimating the vertical ionization energy from the HOMO, 
we present a more direct way to do the calculation. In this calculation, we optimized 
the geometry of both the neutral gas molecule and the positive one charged ion accord¬ 
ingly, for the spin-Unrestricted Hartree-Eock (UHE) wave functions. Then the vertical 
ionization energy is calculated from the ground state energy of positively charged ion 
by subtracting that of the neutral molecule. The results are shown in Table. (|^. 

Eor the electron affinity calculation, we calculated the LUMO with two methods. 
The Hartree-Eock method, computes the closed-shell Restricted Hartree-Eock (RHE) 
wave functions, and the Density Eunctional Theory (DET), computes the closed shell 
densities and Kohn-Sham orbitals in the Generalized Gradient Approximation. The 
results are also shown in Table. ([^. 

The Eig.([^ compares the HOMO (in blue) and LUMO (in green) of R12 on the left 
plot and the HOMO for neutral R12 molecule (in blue) and charged one R12 ion (in 
red) on the right plot. Similar plot for HE01234ze is shown in Eig. As can be seen 
from the HOMO plot, the molecular orbit between the neutral molecule and charged 
molecule are similar, especially for the R12, which has less atoms in the molecule. 
Since the Koopman’s theorem will be exact when the molecular orbit for the neutral 
molecule and charged molecule are the same, we calculated the atom number averaged 
difference between our two methods, which reflect the difference of the molecular 
orbitals. The results are shown in Table. ([8]) 

As the lowest excitation energy can be estimated from the energy gap between the 
HOMO and LUMO, we also make a comparison of these results to those obtained in 
the previous section. Erom the comparison in Table. ([^, we see that, although the ab¬ 
solute values of the excitation energy have relatively large differences between the two 
approaches, the difference of the excitation energies among the molecules considered 
is only marginal. 
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(a) R12 HOMO and LUMO (b) R12 and R12+ HOMO 

Figure 8: A comparison of Rl2{CCl2F2) HOMO (in blue) and LUMO (in green) for 
the neutral molecule (left plot) and the HOMO for neutral R12 molecule (in blue) and 
charged one R12 ion (in red, right plot) . 



(a) HF01234ze HOMO and LUMO (b) HF01234ze and HF01234ze+ HOMO 

Figure 9: A comparison of ll¥01234zt(CFHCHCFs) HOMO (in blue) and LUMO 
(in green) for the neutral molecule (left plot) and the HOMO for neutral HF01234ze 
molecule (in blue) and charged one HF01234ze ion (in red, right plot). 
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7 Estimation of the first Townsend parameter 


Following Ref. l60l . we express the dependance of the first Townsend parameter a as 
a function of the reduced electric field Ej'p 

- = (15) 

P 

where A and B are parameters depending on the gas type and electric field range, and 
p the gas pressure. 

The reference ll60l shows that the first Townsend coefficient at high reduced electric 
field depends almost entirely on the mean free path of the electrons. The mean free 
path, which is defined as 


Am = (16) 

na 

where n is the number of atoms per unit volume and a is the total cross section for 
electron collision with atoms, can be calculated if the environment of gas molecule is 
provided. Firstly, the n is dependent on the gas pressure and can be put by hand. Sec¬ 
ondly, the cf{v), which is a function of electron velocity v, can be calculated from the 
NIST website (SD- Then, noting that at high reduced electric field, the first Townsend 
coefficient can be expressed as a combination of two components: 


(17) 

0(2 =na,E > Em 

where Em is the electric field that drifting an electron over Xm and reach the ion¬ 
ization energy (Jq) of the gas. Therefore, we can estimate with two free parameters A, 
B, which need to be determined from experiments. 

For real gases, we cannot take the free path lengths as a constant and the ionization 
cross section is only a fraction of the total cross section. In that case the estimation 
needs to be modified and use the following equation: 


ai = Ap ’ exp 


a{r) = Ap • exp 



— exp 


Iona \ \ 
eE{r))) 


naiexp 


Iona \ 
eE{r))' 


(18) 


Numerous measurements of the Townsend coefficients are available for standard gas 
mixtures, such as those reviewed in 16^ . 
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Conclusions 


Currently used F-based gases today used in HEP gas detectors are being phased out 
by industry and replaced by eco-friendly substitute gases. This study has reported on 
a general survey of industrially available replacements for HEP gases, discussed their 
physical properties, materials compatibility and safety issues. Parameters of interest 
for their use in HEP detectors have been computed following different approaches 
ranging from parametrizations to quantum chemical calculations: ionization energy, 
electronegativity, number of primary pairs. Statistical methods to compute amplifica¬ 
tion parameters of the ionization shower production such as the Townsend coefficients 
were investigated and preliminary results reported. Promising candidates with lower 
GWP are identified for further studies. 
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8 Appendix - Datasheets 
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Table 9: Gas data sheets (1). 
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Table 10: Gas data sheets (2). 
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Table 11: Gas data sheets (3). 
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